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Abstract. We study both static and dynamic properties of a weakly interacting Bose- 
Einstein condensate (BEC) in a quasi one-dimensional gravito-optical surface trap, where 
the downward pull of gravity is compensated by the exponentially decaying potential of an 
evanescent wave. First, we work out approximate solutions of the Gross-Pitaevskii equation 
for both a small number of atoms using a Gaussian ansatz and a larger number of atoms 
using the Thomas-Fermi limit. Then we confirm the accuracy of these analytical solutions 
by comparing them to numerical results. From there, we numerically analyze how the 
BEC cloud expands non-ballistically, when the confining evanescent laser beam is shut off, 
showing agreement between our theoretical and previous experimental results. Furthermore, 
we analyze how the BEC cloud expands non-ballistically due to gravity after switching off 
the evanescent laser field in the presence of a hard-wall mirror which we model by using a 
blue-detuned far-off-resonant sheet of light. There we find that the BEC shows significant 
self-interference patterns for a large number of atoms, whereas for a small number of atoms, a 
revival of the BEC wave packet with few matter-wave interference patterns is observed. 
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1. Introduction 

Bose-Einstein condensation (BEC) is impossible in a one- or a two-dimensional homogeneous 
system [1, 2], but does occur in atomic traps because the confining potential modifies the 
density of states [3, 4, 5]. Experimentally, a highly elongated quasi-ID regime can be reached 
by tightly confining the atoms in the radial direction, which can be effectively achieved by 
letting the radial frequency be much larger than the axial frequency [6, 7, 8, 9,10, 11,12, 13]. 
However, when the radial length scales become as the order of the atomic interaction length, 
the one-dimensional system can only be described within the Tonks-Girardeau or within 
the super Tonks-Girardeau regime [14, 15, 16], which is experimentally reachable near a 
confinement-induced resonance [17, 18, 19]. 

The interaction of lower dimensional ultra-cold atoms with surfaces has attracted much 
attention in the past few years as their enhanced quantum and thermal fluctuations have 
turned out to play an important role for various technological applications [20, 21, 22, 23]. 
Also under such circumstances, the influence of gravity must be taken into account. Eor 
many years, atomic mirrors were constructed in the presence of a gravitational field by 
using repulsive evanescent waves, which reflect both atomic beams and cold atom clouds 
[24, 25, 26]. The trapping of atoms in a gravitational cavity, consisting of a single horizontal 
concave mirror placed in a gravitational field, is discussed in detail in Refs. [27, 28]. The 
inherent losses of atoms in a gravitational cavity can be reduced by using a higher detuning 
between the evanescent wave and the atomic resonance frequency in a gravitational trap [29]. 
In 1996, Marzlin and Audretsch studied the trapping of three-level atoms in a gravito-optical 
trap by using the trampolining technique without the trampoline [30]. Another approach 
proposed by Saif et al., uses a spatially periodic modulated atomic mirror to yield either a 
localization [31] or a coherent acceleration [32] of material wave packets depending on the 
chosen initial conditions and the respective system parameters. In 2002, Nesvizhevsky et 
al. reported the evidence of gravitational quantum bound states of neutrons in a gravitational 
cavity [33]. Not only is a good confinement geometry necessary for trapping and observing 
the dynamics of atoms, but an experiment also needs an efficient loading scheme loading 
scheme. The experimental group of Rudi Grimm from Innsbruck demonstrated both the 
loading of ^^^Cs atoms [34, 35] and the subsequent creation of a BEC in a quasi-2D gravito- 
optical surface trap (GOST) [36, 37]. More recently, Colombe et al. studied the scheme for 
loading a ®^Rb BEC into a quasi-2D evanescent light trap and observed the diffraction of a 
BEC in the time domain [38, 39]. Later Perrin et at. studied the diffuse reflection of a BEC 
from a rough evanescent wave mirror [40]. 

Motivated by the crucial relevance of gravito-optical surface traps in atomic waveguides 
[41, 42, 43] and atomic chips [44, 45,46, 47], in this paper we study the special case of a quasi 
one-dimensional Bose-Einstein condensate which is trapped orthogonal to the prism surface 
along the vertical axis. In our proposed model the downward pull of gravity is compensated 
by an exponentially decaying evanescent wave (EW), which can be thought of as a mirror 
as it repels the atoms upward against gravity as shown in Eig. 1. In order to deal with the 
hard-wall boundary condition, we apply the mirror solution analogy to the BEC context, and 
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Figure 1. (Color online) Anharmonic GOST potential (solid line) from Eq. (2) in 
dimensionless units, which are explained at the end of the Sec. 2. It consists of a superposition 
of an exponentially decaying optical potential (red circles) due to an evanescent light field 
above a mirror and a linear gravitational potential (dashed-green circles). When the atoms are 
cold enough, they stay in the vicinity of the potential minimum. 


obtain analytical results, which agree with those from numerically solving the underlying one¬ 
dimensional Gross-Pitaevskii equation (IDGPE). Later on, as an interesting applieation, we 
compare our numerical simulation results for a time-of-flight dynamics with the Innsbruck 
experiment for a quasi-2D EEC in a GOST [37]. Surprisingly, our proposed quasi one¬ 
dimensional model agrees even quantitatively with the Innsbruck experiment. Although this 
Innsbruck experiment uses a 2D pancake shaped EEC, when performing the time-of-flight 
expansion vertically the transversely eonfining beam is kept constant, so our quasi-ID model 
for a EEC should apply in this case. 

The paper is organized as follows, the underlying fundamental model for such a quasi- 
ID EEC is reviewed in Sec. 2. Furthermore, we provide estimates for experimentally realistic 
parameters, which we use in our quantitative analysis. Afterwards, we work out approximate 
solutions for the IDGPE wave funetion in the ground state of the system. To this end. See. 
3 performs a modified Gaussian variational ansatz for weak interactions, which corresponds 
to a small number of atoms. For a larger number of atoms, the interaetion strength becomes 
so strong that the Thomas-Fermi solution turns out to be valid, as described in Sec. 4. Then, 
in Sec. 5, we outline our numerical methods and compare them to the previous analytical 
solutions. In Sec. 6, we deal with the time-of-flight expansion of the EEC when the EW is 
removed, showing quantitative agreement with previous experimental results. In Sec. 7 we 
discuss further dynamieal properties of the EEC in a GOST after switching off the evaneseent 
laser field in the presence of the hard-wall mirror. Lastly, we summarize our findings and end 
with brief coneluding remarks. 
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For our ID model of the EEC in a GOST, we assume that we have a dilute Bose gas and that 
the radial frequency is much larger than the axial frequency, i.e. the EEC is cigar-shaped. 
With this assumption, we arrive at the IDGPE [48, 49] 

ih^f{z,t) = + ^(^) ( 1 ) 

On the right-hand side of the equation, the first term represents the kinetic energy of the atoms 
with mass ttib, while the last term describes the two-particle interaction, where its strength 
Gb = 2 N^a^lvjJr is related to the s-wave scattering length aB, and the particle number 
whereas denotes the radial trapping frequency. The anharmonic potential energy V(z) in 
Eq. (1) is produced by both gravity and the exponentially decaying evanescent wave as shown 
inFig. 1[2]: 

V(z) = Voe"'""' -I- mBg;^. (2) 

Here, g is the gravitational acceleration and the constant Vq = rAo/o/(87r^C(53) denotes the 
strength of the evanescent field, where F is the natural linewidth of ^^^Cs atoms, Aq = 852 nm 
is the wavelength of the optical transition, Iq stands for the peak intensity of the EW, and 63 
corresponds to the detuning frequency of the hyperfine sub-level F = 3 of the ^^^Cs atom 
[34, 35, 37]. Furthermore, 1/k = A/2 = XjA.'KsJn? sin^ 6 — 1 represents the decay length, 
where A is the wavelength of the EW, n stands for the refractive index of the medium and 9 is 
the angle of incidence. The potential Eq. (2) has a minimum at z™'’^ = (1/k) ln(VoK/mBg) 
with the axial frequency coz = Note that this potential yields a hard-wall condition with 

< 0) = 00 , because the atoms cannot penetrate the prism, as it is a macroscopic object. 

In order to have a concrete set-up in mind for our analysis, we adapt parameter values 
from the GOST experiments [34, 37]. Eor the EW, we consider the inverse decay length 
to be K = 2/A = 1.43 x 10® m“^, i.e. A ~ 1.4/im. Additionally, we assume an axial 
frequency of 27r x 600 Hz. Eor our atoms in the F = 3 state, the strength of the EW 
is given by Vq ~ 100 x ks pK, where ks is the Boltzmann constant. This potential value 
is within an order of magnitude of the Innsbruck experiments [34, 37]. In view of the quasi- 
ID model, we must satisfy the condition <C cOi, so we assume = 27r x 3 kHz, which 
corresponds to the radial oscillator length k = 0.892 pm. The experiment uses a magnetic 
field for Eeshbach resonance, such that the s-wave scattering length amounts to ob = 440 ao 
with the Bohr radius uq. As both cjz <C cOr and k > are fulfilled, we have, indeed, a quasi 
one-dimensional setting. 

When the atoms are close to the dielectric surface of the prism, we would have to add 
an additional potential contribution due to the van-der-Waals interaction. Here we have to 
distinguish two special cases depending on the distance of the BEC from the surface. The 
regime, where the BEC is close to the dielectric surface, i.e. .2 A/27r, is called the 

Lennard-Jones regime [50, 51]. In the case where the BEC is far away from the surface, i.e. 
2; ^ X/2tt, the regime is called Casimir-Polder regime. In the latter case the Casimir-Polder 
potential can be described within a two-level approximation of the Cesium atoms according to 
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Vcp = —C/^jz^, where the Casimir-Polder coefficient is 6*4 = 1.78 x [52, 53, 54], 

As the EEC does not penetrate very far into the repulsive EW, it is hardly influenced by 
the van-der-Waals potential, which follows from the above values of the GOST experiment. 
Indeed, for the wavelength A = 839 nm [37] and the distance being estimated by the minimal 
distance of the EEC from the surface . 2 ™'" = 4.761 pm., we are in the the Casimir-Polder 
regime. Thus, the value of the Casimir-Polder potential is of the order Vcp = 24.9 x ks pK, 
which is negligible in comparison with the EW potential Vq ~ 100 x ks pK. 
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Figure 2. (Color online) Width A (triangles) and mean position zq (circles) as a function of the 
number of atoms N^. Note that the mean position gives meaningless values for > 3000. 
Thus, the variational ansatz is only successful for quite a small number of particles. 

In view of the forthcoming discussion we use dimensionless parameters as follows. First 
we introduce the dimensionless spatial coordinate 5 = kz . Further, we multiply all terms in 
Eq. (1) by k/( mpg), yielding the dimensionless GPE, 

o El o2 

^) = |“2 ^ ^ (3) 

where the dimensionless kinetic energy constant reads k = (gm|), the dimensionless 

time i = t{mQg) /{hn) and the two-particle dimensionless interaction strength Gb = 
2N^Cjrd^ with Ob = being a dimensionless s-wave scattering length. Additionally, we 
measure energies in units of the gravitational energy m^g/K and get Cjz = hnuJz/gmB as a 
dimensionless frequency, and I^o = n.Vo /gms as a dimensionless strength of the evanescent 
field. 

According to these chosen parameters, the dimensionless quantities have the following 
values. The dimensionless strength of the EW is Vo = 905.7, the dimensionless kinetic energy 
amounts tok = 0.066, the dimensionless s-wave scattering length is given by ob = 0.033, the 
dimensionless radial frequency yields ojr = 1-303, and, finally, the resulting dimensionless 
two-particle interaction strength is Gb = 0.086A^b- From here on, we will drop the tildes for 
simplicity. 
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Figure 3. (Color online) Energy of the EEC as a function of imaginary time for a decreasing 
number of particles from the top to the bottom. 



For the number of particles < 3000, the effective interaction strength is quite small. 
In this limit the EEC has only a reduced extension, so the anharmonic confinement V(z) 
approximately corresponds to a harmonic potential well around its minimum 6.809. 

Therefore, it is reasonable to propose a Gaussian-like ansatz for the wave function in the static 
case. In order to meet the hard-wall condition, however, we modify the Gaussian function 
such that it has the form of a so-called mirror solution [55, 56], 


flz) oc exp 


'-{z - Zof 


'-{z-\-zof' 

2^2 

— exp 

2 A 2 


(4) 


where zq is the mean position and A represents the width. In this way it is guaranteed that 
the wave function meets the hard-wall boundary condition f(0) = 0. In order to find the 
variational parameters zq and A, we minimize the energy of this ansatz following the idea of 
Perez et al. [57, 58]. In order to simplify the expression for the corresponding energy, we 
introduce the parameter 7 = z^jA and normalize the wave function (4) to obtain 


exp(-^jSinh(f) 

f{z) = ^ ’ =. (5) 

y TV^[exp (7^) - 1] 


From this ansatz, we obtain the Gross-Pitaevskii energy 


[Erfc (j - 


2y4e^'7Erf(7) 
2 - 1 ) 



7) + e^^^Erfc (7+7)] 


2 - 1 ) 


- lU 
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2V^A 


( 6 ) 

























Quasi one-dimensional Bose-Einstein condensate in gravito-optical surface trap 


1 


where Erf(y) = e~^^dx = 1 — Erfc (y) denotes the error function. Although, this 

expression cannot be minimized analytically, we can use numerical techniques to extremize 
it with respect to the parameters 7 and A based on the values of k, Gb, and V"o given in Sec. 
2. We can see from Eig. 2 that our variational approach turns out to be valid only for quite a 
small number of atoms. Indeed, the extremization process fails when the condensate has more 
than around 3000 atoms, as then the mean position becomes zero as shown in Eig. 2. Note 
that EEC experiments with such small particle numbers are possible, see for instance Refs. 
[59, 60]. 



Z Z 


Figure 4. (Color online) Probability density plots from (a) Gaussian ansatz (5) and (b) 
numerical calculations for an increasing number of atoms from the top to the bottom. 


4. Thomas-Fermi Limit Solution 

Eor a large enough number of atoms, the effective interaction term and the potential term 
are much larger than the kinetic term. In such a case, an approximate Thomas-Fermi (TF) 
solution is found by neglecting the much smaller kinetic term. Thus, the time-independent 
GPE reduces to 

/i ~ 2 ; + Voe ^ + Gb|'0(2)|^. (7) 

In order to determine the chemical potential, we use the normalization condition, which reads 
in the dimensionless scheme as ^ \ f>f‘ dz = 1. Thus, we get 



where zi and .22 denote the zeros of the integrand. For larger 2; the decaying exponential 
vanishes, and so we have 22 ~ p. Examining the smaller root, we see that for small values 
of 2 and moderate values of p, z/p is quite small, thus a reasonable approximation for this 
root is 2 i fa log (Vq/p). This motivates us to divide the TF solution into two parts: first, the 
soft-wall TF solution, where Vq > p, so that zi is larger than zero. Second, the hard-wall 
TF solution in the case Vq < p, where the lower integration limit zi would be less than zero. 
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so the soft-wall TF wave function would fail. The latter case necessitates to use the mirror 
principle in order to guarantee the hard-wall boundary condition. 


4.1. Soft-wall Thomas-Fermi solution 

With the integration boundaries zi ~ log(Vo//i) and Z 2 ^ p known to be a good 
approximation, we carry out the integration in Eq. (8), yielding 


2 Gb 





(9) 


where we have neglected the small terms {Vq/ p)e~^^ and log {Vq/ p)^ / {2p). Thus, to leading 
order, we have p ~ \/2Gb with a subsequent logarithmic correction. For small changes of the 
chemical potential, the natural log term in Eq. (9) does not change significantly. Therefore 
in Eq. (9) we can substitute \/2Gb for p in the natural log. Using the quadratic equation and 
neglecting small terms in the square root, we obtain the improved approximation 


p ~ \/2 Gb + 1 + log ^Vq/ s/2 Gb^ . 


( 10 ) 


Thus, assuming that Vq > \/2Gb, we obtain the following soft-wall TF solution for 


Zi < Z < Z2. 



( 11 ) 


where p is given by Eq. (10). The function is set to zero for z < zi ^ logiVo/s/^) and 
z > Z 2 ^ p, because the probability density cannot be less than zero. 


4.2. Hard-wall Thomas-Fermi solution 

Eor particle numbers Nq > 2.35 x 10^ the soft-wall TE solution is not valid anymore as 
Zl log(Vo/\/2GB) becomes negative. In order to extend this approximate solution to the 
case where Vq < p, we must work out the corresponding hard-wall TE solution. With the help 
of the mirror analogy [55, 56], we obtain the approximate TE wave function 


f{z) = 


Vo^-, 


41 (1-^-^e 


^ YOp — Z 

Id 11^ 


Gb 7 V iJ- iJ- 


for 0 < z < \zi\ 

for \zi\ < z < Z 2 . 

( 12 ) 


Here, M denotes the normalization constant which is determined from J \f>f dz = 1. Note 
that an analytical derivation of M is not possible, therefore, we performed the respective 
integration numerically. 
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5 . Numerical Methods and Results 

In order to demonstrate the validity of the proposed analytical results, we numerically find the 
ground state of the wave function by propagating the GPE in imaginary time, with the help 
of the split operator technique [61, 62, 63, 64, 65]. For the above mentioned experimental 
parameters and with the value of A^b ranging from 10® to 10® atoms, the ground-state energy 
of the EEC in a GOST does, indeed, quickly converge as shown in Fig. 3. 



Figure 5. (Color online) Comparison of numerical results from GPE (solid lines) with TF 
approximation (circles), (a) A^b = 10^, A^b = 10^, and = 10®. (b) A^b = 10®, A^b = 10^, 
and Nb = 10® from top to bottom, respectively. 

With this technique in mind, we compare our analytical results from Sec. 3 and 4 to 
numerical results and show how the EEC wave function in a GOST changes with increasing 
the number of atoms. For a smaller number of atoms, the variational Gaussian ansatz is 
more suitable as shown in Fig. 4, whereas for a larger number of atoms the Thomas-Fermi 
approximation turns out to be quite accurate as shown in Fig. 5. The variational Gaussian 
ansatz roughly reproduces the mean, but it is rather symmetrical, unlike the numerical 
results as the number of atoms approach 2400, see Fig. 4. Qualitatively, the EEC width is 
proportional to the number of atoms in a GOST. However, due to the EW decaying exponential 
potential, the EEC cannot expand in the negative 2 :-direction, so the EEC wave function takes 
up a triangular shape for A^b larger than 10® as shown in Fig. 5. The agreement between 
numerical and analytical TF results is much better for larger number of atoms. Note that the 
EEC wave function is quite wide for 10® atoms in Fig. 5. 

The van-der-Waal forces with the surface can demolish the EEC, so it is necessary to 
have a larger EW potential. Therefore, keeping in mind recent developments in the laser field 
technology, it is possible to increase the EW strength Vq by increasing the laser power [66]. 
Thus, we explore now the parameter space of our model. Irrespective of the number of atoms, 
the mean position of the EEC in GOST increases only moderately with Vq. For the soft-wall 
case Vq > p, as shown in Fig. 6, the maximum of the wave function occurs at log(Vo), but for 
the hard-wall case Vq < p the maximum of the EEC wave function exists at |log(V[)/\/2G^) |. 
Thus, due to the interaction term, for a large number of atoms, the maximum of the wave 
function no longer remains within the minimum of the trapping potential ~ log(Vo). 
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Figure 6. (Color online) Mean position of the BEC versus EW strength Vq for decreasing 
number of atoms from the top to the bottom. 


The BEC wave function in a GOST becomes asymmetric for larger interaction 
strengths. Therefore we quantify the BEC width based on the standard deviation a = 
\/< > — < z >^, where < • >= f • \'ip{z)f dz denotes the expectation value. As shown 
in Fig. 7, the BEC standard deviation grows extremely rapidly with increasing number of 
atoms A^b- However, on the other hand, changing Vq only slightly affects the standard 
deviation a. 
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Figure 7. (Color online) Standard deviation a of BEC wave function increases with inter¬ 
particle interaction strength Gb for EW strength Vq = 906. 


6. Time-of-flight (TOF) Expansion 

The standard observation of a BEC is based on suddenly turning off the trapping potential 
and allowing the atoms to expand non-ballistically. The resulting time-of-flight (TOF) 
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measurements are performed either by acquiring the absorption signal of the probe laser beam 
through the falling and expanding BEC cloud, or by measuring the fluorescence of the atoms 
which are excited by a resonant probe light [40]. 



t (ms) 


Figure 8. (Color online) Fraction of remaining atoms during time-of-flight during a vertical 
expansion for Vq = 453. Initially, total number of atoms is Nb = 2400. Here the circles stem 
from the Innsbruck experimental [37], whereas the solid line shows our numerical results. 

In the Innsbruck experiment, the remaining number of atoms is measured after allowing 
atoms from the GOST to expand vertically by suddenly turning off the EW [37]. Note that 
some particles are lost due to thermalization processes which occur when the particles hit the 
prism or due to the van-der-Waal forces with the surface. Although this Innsbruck experiment 
uses a 2D pancake shaped BEC, when performing this vertical expansion the transversally 
confining beam is kept constant, so our quasi-ID model for BEC should apply in this case. 
Using the experimental parameters in Ref. [37], we numerically reproduce their vertical 
expansion curve (their Figure 2), as shown in Fig. 8. To this end, we use A^b = 2400 atoms and 
Vo = 453, yielding an initial condensate wave function with dimensionless standard deviation 
a = 0.86 and the dimensionless mean position zq = 6.40. We propagate this wave function 
without the hard-wall boundary condition. Then we approximate the fraction of remaining 
atoms by t)\'^dz, as atoms in the BEC wave function extending past z = 0 are lost by 

sticking to the surface. As the interaction term is quite small in the Innsbruck experiment, the 
standard deviation of the BEC is small and remains roughly constant during the time-of-flight, 
so the loss of atoms does not affect time-of-flight expansion significantly for t < 0.7 ms as 
shown in Fig. 8. 

We also simulated the TOE expansion without the prism at z = 0 for different particle 
numbers A^b- We see that the mean position of the BEC drops due to gravity as shown in Fig. 9 
with different rates, which strongly depend upon the number of confined atoms. At the same 
time the BEC width, which is proportional to the standard deviation, increases according to 
Fig. 9. 
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Figure 9. (Color online) (a) Mean position zq and (b) standard deviation a of BEC density 
in the time-of-flight expansion for decreasing number of atoms from the top to the bottom for 
Vo = 906. 


7. Dynamics of BEC on hard-wall mirror 

Concerning the dynamics of the BEC in a GOST, we consider two cases. First we describe the 
dynamies of the BEC with A^b = 1800 after switehing off the evanescent potential and letting 
the BEC atoms fall on a hard-wall mirror. The latter could be experimentally realized by a 
blue-detuned far-off-resonant sheet of light, and is modeled theoretically via the boundary 
condition E(0) —> oo at f > 0. Thus, the BEC of A^b = 1800 atoms has the dimensionless 
mean position zq = 10.7 and the dimensionless standard deviation a = 4.28 at time t = 0. 
We observe the matter-wave interference pattern formed upon releasing the eondensate from 
GOST as shown in Fig. 10, as atoms impinging on the ‘hard-walT at the origin ‘bounce’ 
baek. For short times the atoms remain near the hard-wall surface, so the BEC dynamies is 
characterized by the reflection of atoms from the hard-wall mirror. But for larger times those 
atoms, whieh are far away from the hard-wall mirror, are reflected above the hard-wall mirror 
due to collisions among themselves as shown in Fig. 10. The number of atoms A^b = 1800 
is so large that the BEC is staying quite close to the mirror, so we have not seen any total 
reflection of the BEC wave packet. 

In order to see the latter scenario, we need a small number of A^b = 30 particles far away 
from the hard-wall mirror. This is realized by the EW potential whieh could be 

implemented by trapping atoms in a MOT above the surface, so that, once the EW trap is 
switehed off, the atoms have enough momenta when they hit the hard-wall mirror as shown in 
Fig. 11 . Similar to experiments of photonic bouncing balls [67] and plasmonic paddle balls 
[68], the BEC shows significant self-interference patterns, for example in the time intervals 
t=24-34 and 83-91 as shown in Fig. 11 , whieh originate from the interferenee of ineoming and 
reflecting BEC wave packets. It is worth mentioning that a smaller initial width of the BEC 
wave paeket would lead to finer revivals and a larger initial width of the BEC would lead to 
larger interference regions. The evolution of a BEC falling under gravity and bouncing off a 
hard-wall mirror formed by a far-detuned sheet of light was already observed experimentally 
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Figure 10. (Color online) Dynamics of BEC in presence of gravity with hard wall at 2 : = 0 
for A^b = 1800 number of atoms. Here the color scale represents the density of the BEC. 


by Bongs et al. in both the soft-wall and the hard-wall regime [69]. In the soft-wall regime, 
they have recorded that the BEC is reflected up to three times off the optical mirror in the 
lossy environment. Due to a large two-particle coupling strength, which in turn results in a 
condensate with a larger width, this group also observed a splitting of the BEC into two parts 
close to the upper turning point of the BEC. This effect is heuristically modeled by a GPE 
dynamics with assuming that the two-particle interaction strength decreases exponentially in 
time. In our case we restricted ourselves to the evolution of the BEC with a constant coupling 
constant, so we did not observe any splitting of the BEC in our simulation, but we do observe 
the BEC for longer time intervals as shown in Fig. 11. In our simulation, we observed more 
than three reflections of the atomic cloud in a lossless environment. The center-of-mass of 
the BEC wave packet shows periodic revivals with the dimensionless time period t 58.2 as 
shown in Fig. 11 . Quantitatively, the classical revival time can be determined in dimensionless 
units as f = 2 ^/2z^Jk = 56.9, where zq represents the mean position of the BEC and the 
dimensionless energy constant k is defined below Eq. (3). 

8. Conclusion 

In summary, we studied the behavior of a Cs BEC in a quasi-ID optical surface trap. We 
have developed approximate solutions to the GP equation for both small and large numbers 
of atoms. In the former, we have used the variational ansatz technique, while in the latter we 
have used the Thomas-Fermi approximation. Later on, we compared the analytical approach 
with numerical results, which agreed quite well for a wide range of atom numbers A^b- 

Furthermore, we have numerically reproduced the experimental result of Ref. [37], 
where a 2D BEC is confined in the radial direction, but is allowed to expand in the vertical 
direction freely. This indicates that our analysis could be extended beyond the ID case. Our 
model suggests that for a small particle number A^b the BEC retains its Gaussian shape in the 
expansion and falls due to gravity. As suggested by Fig 9, for larger number of atoms, the 











Quasi one-dimensional Bose-Einstein condensate in gravito-optical surface trap 


14 


tsj 



Figure 11. (Color online) Numerical results for the BEC density with initial 

dimensionless mean position zq = 27.2 and dimensionless standard deviation a = 0.99 for 
A^b = 30 number of atoms. BEC experiences full periodic revivals, however incoming and 
reflecting BEC wave packets lead to larger matter-wave interference regions at time intervals 
t=24-34 and 83-91. Here the color scale represents the density of the BEC. 


Standard deviation does not expand as fast as compared to small number of particles, therefore 
we conclude that the initial number of the particles plays a significant role in the expansion 
of the BEC cloud. Afterwards, we investigated the dynamics of the BEC in the presence of 
gravity and a hard-wall boundary condition, where we observed self-interferences and revivals 
of the wave packet. The observation of the bouncing of the BEC can be used to characterize 
and determine mirror properties such as roughness and steepness. All our results can be 
applied to develop atomic interferometers for a BEC. 
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